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ABSTRACT 

That the stellar halo of the Milky Way has a density profile which to first approxima- 
tion satisfies p oc has been known for a long time. More recently, it has become clear 
that M31 also has such an extended stellar halo, which approximately follows the same 
radial scaling. Studies of distant galaxies have revealed the same phenomenology. Also, 
we now know that the density profiles of the globular cluster systems of our Galaxy and 
Andromeda to first approximation follow p oc r^^, E oc R^^ in projection. Recently, 
diffuse populations of stars have been detected spherically surrounding a number of 
Galactic globular clusters, extending much beyond the Newtonian tidal radii, often 
without showing any evidence of tidal features. Within the standard Newtonian and 
GR scenario, numerous and diverse particular explanations have been suggested, in- 
dividually tailored to each of the different classes of systems described above. Here 
we show that in a MONDian gravity scenario, any tenuous halo of tracer particles 
forming a small perturbation surrounding a spherically symmetric mass distribution, 
will have as an equilibrium configuration precisely a p oc r^^ scaling. 
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1 INTRODUCTION 

Progress in the debate between the hypothetical physical 
reality of astrophysical dark matter and the option of modi- 
fying gravity at the low acceleration regime, will hinge upon 
the exploration of as many independent lines of enquiry as 
possible. Whereas the rotation curves of galaxies can be 
adequately reproduced from either point of view (Milgrom 
1983, Milgrom 1994, Sanders & McGaugh 2002, Swaters et 
al. 2010), a variety of recent studies have shown results in 
tension with the standard scale-invariant gravity plus dark 
matter paradigm, and being more in line with generic MON- 
Dian gravity approaches. 

We shall use the term MONDian to refer to modi- 
fied gravity theories in which, in the low velocity limit, the 
force per unit mass between a test particle and a spherical 
mass distribution will shift from the Newtonian expression 
of GM/r^ to an {GMaoY^'^ /r behaviour for acceleration 
scales below Milgrom's ao, independently of the fundamen- 
tal underlying theory of gravity which might lead to such 
a behaviour e.g. Bekenstein (2004), Moffat & Toth (2008), 
Zhao & Famaey (2010), Bernal et al. (2011), Mendoza et al. 
(2011), Capozziello & De Laurentis (2011). The most salient 
features of such schemes are equilibrium velocities which be- 
come independent of distances at a value ~ (GA/ao)^''^. 



In Lee & Komatsu (2010) and Thompson & Nagamine 
(2012) it has been shown that the infall velocity of the two 
components of the Bullet cluster is incompatible with expec- 
tations of full ACDM predictions, and is actually a challenge 
to the standard gravity theory, as it surpasses the escape ve- 
locity of the combined system. Recently, Kroupa (2012) has 
shown that tidal dwarf galaxies, which under the standard 
scenario are transient tidal clumps of galactic material hav- 
ing out of equilibrium dynamics, actually fall on the same 
TuUy-Fisher relation as all dwarf galaxies. This appears as 
an uncanny coincidence from the standard gravity perspec- 
tive, where the dynamics of normal dwarfs are thought to 
be determined by their dominant dark matter halos. The 
result is expected under MONDian gravity schemes, where 
once the baryonic mass of the system is fixed, the dynamics 
will uniquely follow, as it is indeed observed. 

Along the same lines, in Scarpa et al. (2011), Hernandez 
& Jimenez (2012) and Hernandez et al. (2012b) it has been 
shown that the velocity dispersion profiles of a number of 
Galactic globular clusters stop falling radially along Keple- 
rian expectations and settle to finite asymptotic values on 
crossing the a < ao threshold. The standard gravity expla- 
nation for these profiles, that it is the tides of the Milky Way 
that dynamically heat the outskirts of the clusters observed 
e.g. Lane et al. (2010), appears suspect, as the Newtonian 
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tidal radii can be shown to exceed the points where the ve- 
locity dispersion profiles flatten by large factors, and because 
the total globular cluster masses and asymptotic velocity 
dispersion values follow the galactic TuUy-Fisher relation, as 
expected under MONDian gravity schemes. Finally, in Her- 
nandez et al. (2012a) we showed that the relative velocities 
of extremely wide binaries do not follow the expectations of 
full galactic dynamical Newtonian simulations, but diverge 
from the Keplerian decline with radius to settle also at fi- 
nite relative velocities on crossing the a < ao threshold of 
MONDian gravity proposals. 

In this letter we show that the density profile of a popu- 
lation of trace particles surrounding a spherical mass distri- 
bution in MONDian gravity will naturally follow a p oc 
profile. Under standard gravity approaches, the ubiquitous 
nature of p oc profiles has to be addressed through a 
variety of highly specific explanations, individually tailored 
to each of the diverse classes of systems where these profiles 
have been observed. Examples of the above are mergers and 
tidal dissolution of accreted sub-structure for the tenuous 
stellar halos surrounding our Galaxy, M31 and also the re- 
cently detected ones around external galaxies (e.g. Bullock 
et al. 2001, Abadi 2006), and the compression of tidal tails 
at apocentre or disk shock heating for the " extra tidal" fea- 
tures surrounding Galactic globular clusters, e.g. Da Gosta 
(2012). 

The above situation contrasts with the appearance of a 
direct equilibrium solution for p oc profiles under MON- 
Dian gravity schemes which we show here, and which we 
suggest is one further piece of evidence pointing in the di- 
rection of the necessity of a modified gravity regime at low 
acceleration scales. 

In Section (2) we derive the first order predictions for 
equilibrium tenuous stellar halos under a MONDian gravity 
scheme, and also under Newtonian gravity within dark mat- 
ter halos, and under Newtonian gravity in the absence of 
dark matter halos. In section (3) we compare the theoretical 
estimates with the observed profiles for a variety of systems, 
finding the MONDian prediction a good fit to the observed 
situation, over a wide variety of scales and classes of astro- 
nomical systems. Section (4) presents our conclusions. 



2 FIRST ORDER DYNAMICAL 
EXPECTATIONS 

We shall model the physical situation which applies generi- 
cally for a MONDian scenario, where in the a « ao limit, 

the force between a test particle and a spherically symmet- 
ric mass distribution becomes — (GM(r)ao)^^^/r, when one 
introduces no modifications to the Newton's second law e.g. 
Hernandez et al. (2010), Mendoza et al. (2011). 

Assuming spherical symmetry, and taking the deriva- 
tive of the kinematic pressure, the equation of hydrostatic 
equilibrium for a polytropic equation of state P = Kp'^ is: 



djKp'' 
dr 



In going to isothermal conditions, 7 = 
Binney & Tremaine (1987), and we get 

— — = -V0. 
p dr 



(1) 

1 and K = (T^. e.g. 

(2) 



writing p 
written as: 



(47rr^) ifiMIll the above equation can be 



/ dM{r) 
\ dr 



d^M(r) 
dr^ 



= -Wct>{r) 



(3) 



where a is the isotropic Maxwellian velocity dispersion for 

the population of stars. The above treatment is common, 
and can be found in e.g. Hernandez et al. (2010), where we 
used it in the modeling of dSph galaxies, systems charac- 
terised by flat velocity dispersion profiles, obtaining mass 
models consistent with observed velocity dispersion, half 
mass radii and total masses, in the absence of dark mat- 
ter. Another example of a similar treatment is found in Her- 
nandez & Jimenez (2012), where we modelled eight GGs in 
our Galaxy, accurately reproducing all observational con- 
straints. 

As an illustrative example we can take V(A('r) = 
GM{r)/r'^ for the right hand side of equation (3), the Newto- 
nian expression appearing for a >> ao. Looking for a power 
law solution for M(r) = Mq{t / ro)'^ , we get: 



GMo 



(4) 



and hence m = 1, the standard isothermal halo, M(r) = 
2(T^r/G, having a constant centrifugal equilibrium velocity 
= 2a^ and infinite extent. In going to the MONDian limit 
of a « ao, V0(r) = {GM{r)aoY''^ /r, equation (3) yields: 



m — 3 



n/2 



(5) 



In this limit m = 0, we obtain M(r) = Mo and = 3cr^ = 
[GMoaoY^'^ 1 the expected Tully-Fisher scaling of the cir- 
cular equilibrium velocity with the fourth root of the mass, 
with rotation velocities which remain fiat oven after the mass 
distribution has converged, thus, rigorously isothermal halos 
are naturally limited in extent. It is interesting that in this 
limit the scaling between the circular rotation velocity and 
the velocity dispersion is only slightly modified as compared 
to the Newtonian case, with the proportionality constant 
changing from 2 to 3, for the squares of the velocities. 

We can now look for the behaviour of a tenuous stel- 
lar halo in the MONDian regime, and hence at large dis- 
tances, around a mass distribution which has essentially 
converged, by looking at eq.(2) and writing the R.H.S. as 
— {GMtotaoY^'^/r, where Mtot is the total mass of the galac- 
tic or stellar system. The p in the L.H.S. of this same equa- 
tion now refers to the density distribution of essentially test 
particles making up a trace population, e.g. a stellar galac- 
tic halo, the globular cluster distribution around a large 
galaxy, or the faint halos of "extra-tidal" stars surround- 
ing Galactic globular clusters. Using also the result of eq.(5) 
of 3cr^ = {GMoaof^, eq.(2) yields: 

^ = -3^, (6) 
dr r 

which can then be integrated directly to yield: 

p(r) = po{ro/rf. (7) 

Such a volumetric distribution can be projected anar 
lytically along one direction to yield the projected surface 
density profile 
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S(i?) = 



2E? ■ 



(8) 



Of course, p oc is only an approximation which will 
be valid over a limited radial range. In fact, mass profiles 
for isothermal solutions converge to finite total masses and 
radii, as shown in e.g. Hernandez ct al. (2010), Hernandez 
& jimcncz (2012). In closer detail, the density profiles will 
steepen as one moves further out e.g. as seen in the broken 
power law fits for the Galactic stellar halo reported by Sesar 
et al. (2011). 

Under the Newtonian expression of V(t> = GM{r) /r^ , 
the equivalent development for a trace population in the 
halo of a galaxy having the same rotation curve as the one 
leading to eq.(6), M(r) = 2cr^r/G, where this time M{r) 
refers mostly to the hypothetical dark matter component, 
yields: 



p(r) = po{ro/r)'^ 



(9) 



as the expression corresponding to cq.(7). The case of a 
trace population around an essentially converged total mass 
in Newtonian dynamics, e.g. the tenuous stellar halos sur- 
rounding the Galactic globular clusters, V0 
yields: 



GMtot/r^ 



p{r 



poe 



(GMtot/o-^r) 



(10) 



a density distribution which tends to a constant at large 
radii. Thus, we sec that equilibrium configurations of 
isothermal tracer populations in the MONDian regime will 
follow p{r) oc density profiles, while under Newtonian 
gravity the same populations within the corresponding dark 
matter halos will show much shallower isothermal profiles, 
which in the absence of dark matter halos e.g. tenuous stel- 
lar halos about globular clusters, will have density profiles 
as given by eq.(lO). 



3 OBSERVATIONAL COMPARISONS 

In this section we review the observational situation of ten- 
uous tracer population halos, which now spans a very wide 
range of systems and astrophysical scales. We begin with a 
number of recent determinations of the density structure of 
the stellar halo of our Galaxy. Morrison et al. (2000) im- 
plement a careful disk/halo star separation criteria, and ob- 
tain p(r) oc for the stellar halo of the Milky Way. Also, 
Battaglia et al. (2006) obtain an almost flat velocity disper- 
sion profile for the halo stars in the Milky Way, providing 
empirical support for the use of constant a profiles in the 
theoretical description. Juric et al. (2008) report a single 
power law fit p(r) oc r~^-^^°-^, while Bell et al. (2008) find 
halo profiles having more structure than simple power laws 
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Suprime-Cam instrument on the Subaru telescope found 
T,{R) oc j?-2 "±o is^ and lastly Gilbert et al. (2012) mea- 
sure S(_R) oc 7^-2. 2±o. 2 very large distances, 175 kpc, 
coming to T,{R) oc R-^ o±o.5 20fcpc < R < 90kpc once 
the kiematical substructure is removed . It appears that the 
stellar halo of M31 is a classic example of the tenuous pop- 
ulations described by eq. (7). 

Regarding stellar halos of more external galaxies, tenu- 
ous extended stellar halos have been detected over the past 
few years surrounding numerous systems. Recent detections 
include Cockcroft et al. (2012) who find evidence for an ex- 
tended stellar halo about M33 and Bakos & Trujillo (2012) 
who report finding such structures about a sample of 7 late- 
type spirals from the SDSS, in all cases with total masses 
amounting to only a few percent of the total baryonic mass 
of the host galaxies. Although the above observations can 
not yet yield secure projected density profiles, Jablonka et 
al. (2010) report a tenuous stellar halo about NGC 3957 with 
a projected scaling E(i?) oc 7^-2.76±o.43^ while Baker et al. 
(2009) find a faint stellar halo about M81 with a projected 
scafing S(i?) oc and Bailin et al. (2011) observe a 



stellar halo about NGC 253 with E(i?) oc i?" 



In going 



to larger samples, Zibetti et al. (2004) showed through the 
stacking of images from 1047 edge-on spiral galaxies from 
the SDSS, that these very generally present extended tenu- 
ous stellar halos with volumetric radial density profiles well 
described by p(r) oc r~^. 

In going to the spatial distribution of a different tracer 
population, this time the globular cluster systems of galax- 
ies, it has been well know for many years that the density 
profile of the GC system of the Milky Way very accurately 
follows a p(r) oc profile, e.g. Racine & Harris (1998), 
Surdin (1994). Looking in more detail, more recently Bica 
et al. (2006) find a p(r) oc r"" profile with 3.2 < n < 3.9 for 
all Galactic globular clusters, while the metal rich popula- 
tion follows also a power law, this time with p(r) oc .f^'^-2±o.2 
for large radii, and p(r) cx f-3-2±o.9 j£ includes the effects 
of oblateness in the distribution. The globular cluster sys- 
tem of M31 has a projected power law density profile also in 
agreement with the expectations of eq.(7), e.g. Racine (1991) 
determined S(7i) oc R~^. More recently and in more detail, 
Huxor et al. (2011) obtained a best fit profile composed of 
three distinct power laws, which however, if modelled as a 
single power law for r > Ikpc, can be approximated by the 
same S(i?) oc R~^ law found earlier by Racine (1991). 

Going to more external galaxies, Perelmuter & Racine 
(1995) found a best fit E(7?) oc R'^ scafing for the glob- 
ular cluster system of M81. Harris et al. (1984) found the 
outer projected radial distribution of globular clusters in 
NGC 4594, the Sombrero galaxy to be well described by a 
E(i?) oc R~^ profile. Harris and van den Bergh (1981) found 



to yield better fits, but still, p(r) oc r for the preferred also p(r) oc r scalings for the globular cluster systems 



single power law model. Finally, Sesar et al. (2011) obtain 
a broken power law as the most accurate description, but 
again, a best fit single power law of p(r) oc r~^'^. There is 
clearly a broad radial range over which the best fit single 
power law model yields a slope as expected from eq. (7). 

Recent detailed studies of the stellar halo of Andromeda 
using a variety of techniques and data from the largest 
modern facilities have reached a consensus for a p oc 
structure. Ibata et al. (2007) obtain T,{R) oc i?-i-9i±o ii 
with data from the CFHT, Tanaka et al. (2010) using the 



around 7 elliptical galaxies. More recent studies have found 
a spread in the power law slopes of projected density profiles 
for globular cluster systems, but taken as a whole, "typical 
projected power- law indices range from 2 to 2.5 for some 
low-luminosity Es to 1.5 or a bit lower for the most massive 
giant ellipticals." to quote from the review from Brodie & 
Strader (2006). 

Since the studies of Grillmair et al. (1995) and Leon et 
al. (2000), a number of tenuous stellar halos associated to 
Galactic globular clusters have been detected. The problem 
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of determining structural parameters is harder than in the 
cases of the stellar halos surrounding galaxies, as the overall 
numbers of stars arc much lower, and the problem of con- 
tamination by foreground and background sources, as well 
as by obscuration, is significant. More modern studies have 
found a large range of power law slopes in the outskirts of 
globular clusters, McLaughlin & van der Marel (2005) find 
projected indexes going from -2 to -6, but do report that a 
population of the most mgissive clusters shows indexes close 
to oc R~^, and conclude that the extended halos en- 

veloping the clusters they study are suggestive of a generic 
equilibrium feature, rather than being transient structures. 
Jordi & Grebel (2010) report projected power law indexes 
for tenuous stellar halos surrounding 17 Galactic globular 
clusters; their most reliable results span values from -1 to 
-4. These authors also note features which are problematic 
for a standard gravity interpretation, in the cases of e.g. 
NGC 7089 and Pal 1, whose stellar halos arc clearly spheri- 
cal with no sign of any tidal features, in spite of extending 
in both cases much beyond their Newtonian Jacobi radii. 
Carballo-Bello et al. (2012) perform a similar study (also in- 
cluding careful CMD modelling to limit contamination) for 
the extra-tidal halos of 19 Galactic globular cluster. Again, 
reported projected power law indexes span a broad range 
from close to -2 to about -4, excluding clusters showing clear 
tidal features. Note that for two of the three clusters which 
overlap with the Jordi & Grebel (2010) sample, NGC 4147 
and NGC 5272, Carballo-Bello et al. (2012) report power law 
indexes of -2.8t'^ °l and -3.18i^:^g, while Jordi & Grebel 
(2010) assign to these same clusters values of —1.48 ± 0.24, 
and of —0.94 ± 0.38 respectively, for a comparable radial 
range. This simply illustrates that the observational situa- 
tion is far from converging to definitive answers regarding 
these systems. 

Further, in the case of globular clusters the problem is 
intrinsically less clear than for the galactic stellar halos, as 
internal dynamical evolutionary effects might play a part, as 
well as the gravitational perturbations due to the crossing 
of the Galactic disk. It is however clear that Galactic globu- 
lar clusters often, if not always, are surrounded by tenuous 
stellar halos, which from the Newtonian point of view, must 
be regarded as "extra tidal" structures. The smooth and 
round appearances often observed, with all absence of tidal 
tails, hence become a problem. This problem does not ap- 
pear under MONDian gravity, where satellite systems are 
generically expected to be much more robust to tides, e.g. 
Hernandez & Jimenez (2012). Interestingly, Mackey et al. 
(2010) find a E(i?) oc power law structure for the tenu- 
ous stellar halo surrounding an extremely isolated globular 
cluster in M31, which begins as n « 2.5, and then breaks 
further outwards to n fs 3.5. Also, note that most of the 
S(i?) oc R~" indexes in the two recent Jordi & Grebel (2010) 
and Carballo-Bello et al. (2012) studies cluster about n = 3. 
The latter authors show that the group of clusters showing 
clear tidal features have profiles much steeper than those 
lacking them. The presence of as yet undiscovered tidal fea- 
tures among these clusters would tend to artificially steepen 
their profiles. Note also that already Grillmair et al. (1995) 
cautioned that the difficulties of background subtraction and 
obscuration corrections will lead to systematics which tend 
to yield overestimates in n. On the other hand, the large 
indexes sometimes reported for globular clusters could also 



be detections of the steepening in the profile expected under 
MONDian gravity models on approaching the final radius, 
already mentioned in the discussion following cq.(8). 

It thus appears clear that extended tracer population 
halos, stars or globular clusters, having a small fraction of 
the light of their host systems, galaxies or globular clusters, 
are a common feature. Also, such halos are generally never 
far from the predictions of eq.(7) for equilibrium configu- 
rations of isothermal tracer populations in the MONDian 
regime, precisely p(r) oc or T.{R) oc R-'^. Within the 
standard gravity interpretation, explanations to the power 
law density profiles of galactic stellar halos have been pro- 
posed in terms of the accretion and tidal dissolution of sub- 
structure falling into the main galaxy e.g. Bullock et al. 
(2001), Bullock & Johnston (2005), and Abadi et al. (2006). 
However, the overall smoothness and uniformity in stellar 
properties of these systems has been pointed out as problem- 
atic for the standard explanation, which naturally implies a 
degree of randomness in the accreted material, e.g. Ibata et 
al. (2007), and Bell et al. (2008), who also find from simula- 
tions stellar halos not matching observations in terms of the 
substructure details. Ibata et al. (2007) also show that stan- 
dard simulations sometimes yield exponents of p{r) oc r~" 
inconsistent with observations, with n = 4 or even n = 5. 
Also, a further explanation must be sought for the observed 
profiles of the globular cluster systems surrounding galax- 
ies, and yet another for the remarkably smooth "extra-tidal" 
stellar halos surrounding many of the globular clusters of the 
Milky Way. 



4 CONCLUSIONS 

We have shown that under a MONDian gravity force law, the 
density profiles of tenuous tracer population halos surround- 
ing spherical mass distributions will be well approximated 
by p oc scalings. 

We suggest that such equilibrium configurations provide 
a more natural, and certainly more general, explanation for 
the observed close to p oc r^"^ behaviour of; the stellar halos 
surrounding the Milky Way, M31, and a variety of external 
galaxies, the density profiles of the globular cluster systems 
in our Galaxy and Andromeda, and the radial structure of 
the "extra tidal" stellar halos recently observed surrounding 
a number of Galactic globular clusters. 
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